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PREFACE 


This  report  documents  the  aerodynamic  coefficients  and  stability  derivatives 
resulting  from  a  series  of  free  flight  tests  of  a  generic  missile  configuration  with  circular 
and  non-circular  cross  sections  over  a  Mach  number  range  of  0.75  through  3.5.  All 
configurations  had  a  common  rectangular  tail  fin  design  and  were  either  of  a  three-fin  of 
four-fin  variety  depending  upon  the  missile  cross-section  These  tests  were  conducted  in 
the  USAF  Aeroballistic  Research  Facility,  located  at  Eglin  AFB,  FL.  The  period  of 
testing  covered  a  time  period  of  1997  to  2000. 
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Research  Laboratory  Munitions  Directorate,  Eglin  Air  Force  Base,  Florida  32542-5434. 
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conducted  the  test  that  included  launch,  instrumentation,  data  acquisition,  and  image 
processing. 
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SECTION  I 


INTRODUCTION 


Current  munition  designs  are  departing  from  conventional  right  circular  cross- 
sections  for  many  reasons  including  enhanced  performance,  increased  range,  increased 
maneuverability,  and  stealth.  Improved  manufacturing  and  production  capability  for  both 
metal  and  composite  structures  is  also  leading  weapon  designers  to  consider  non-circular 
vehicle  shapes.  The  objective  of  this  effort  is  to  obtain  an  experimental  aerodynamic 
database  on  several  basic  cross  sectional  shape  variations  within  a  free-flight  ballistics 
range.  The  results  will  serve  to  supplement  existing  data  and  will  be  used  to  improve 
current  analytical  design  methodologies  that  estimate  vehicle  aerodynamics.  A  total  of 
seven  configurations  were  tested.  The  models  were  designed  to  have  an  equivalent  cross 
sectional  area  and  identical  fin  planform  to  isolate  effects  due  to  the  body  alone.  The 
primary  goals  for  this  effort  are  (1)  to  determine  shape  effects  on  the  aerodynamics  and 
(2)  provide  an  experimental  database  for  aeroprediction  methodology. 

The  free-flight  ballistic  range  has  several  advantages  compared  to  wind  tunnels 
for  aerodynamic  research.  The  most  important  one  is  that  the  test  object  is  in 
unrestrained  flight.  Thus,  no  model  support  (sting)  or  wall  interference  effects  are  present 
during  the  measurement  of  the  data.  Additionally,  the  ballistic  range  allows  for  the 
determination  of  aerodynamic  stability  coefficients  and  derivatives  that  are  not  easily 
measured  in  a  wind  tunnel.  Each  trial  conducted  in  the  ballistic  range  results  in  a  unique 
initial  starting  attitude  (initial  conditions)  upon  launch  resulting  in  pitch  and  yaw  motion 
that  enables  the  test  engineers  to  determine  the  dynamic  derivatives  and  coefficients. 
Analyzing  two  or  more  trials  simultaneously  results  in  a  common  set  of  aerodynamic 
parameters  of  the  configuration  independent  of  initial  conditions.  This  can  be  done  over 
a  range  of  Mach  numbers  to  determine  the  aerodynamic  parameters  as  a  function  of  Mach 
number.  It  is  also  desirable  to  have  trials  with  sufficient  pitch  and  yaw  amplitudes  to 
cover  the  angle-of-attack  regime  of  interest.  A  drawback  of  ballistic  range  testing  is  size 
limitations.  Typical  free-flight  models  are  sub-scale  so  matching  the  full  scale  Reynolds 
number  if  difficult  or  impossible.  However,  this  is  also  a  limitation  of  wind  tunnel  tests 
and,  in  either  case,  care  must  be  used  when  interrupting  results. 

One  challenge  in  free  flight  testing  is  the  model  design  and  construction.  Since 
the  model  is  typically  launched  via  circular  tube  (e.g.,  powder  gun)  it  must  be  encased 
within  a  sabot  for  launch.  In  addition,  the  model/sabot  package  must  be  sized  to  fit  the 
particular  launcher  and  be  capable  of  withstanding  the  in-bore  accelerations.  Therefore, 
the  size  of  the  model  is  often  small  and  results  in  high  precision  tolerances  during 
manufacture.  The  in-bore  acceleration  requirement  means  that  the  model/sabot  design 
must  capable  of  surviving  the  launch  cycle.  Hence,  the  models  require  adequate  strength 
to  be  launched  without  suffering  structural  damage. 
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This  report  documents  the  aerodynamic  coefficients  and  stability  derivatives 
extracted  from  trajectory  data  collected  during  free  flight  tests  of  five  non-axisymmetric 
configurations.  In  addition,  two  circular  configurations  were  tested  to  establish  a  baseline 
reference  for  the  results.  The  models  were  designed  so  that  the  cross-sectional  area  was 
common  for  all  configurations.  The  tests  were  conducted  in  the  USAF  Aeroballistic 
Research  Facility.  During  this  period  of  the  tests,  the  Flight  Vehicle  Integration  Branch 
of  the  Air  Force  Research  Laboratory  Munitions  Directorate  solely  maintained  and 
operated  the  Facility.  Recently,  the  University  of  Florida  Graduate  Engineering  and 
Research  Center  (UF/GERC)  has  become  a  partner  in  the  management  of  the  facility. 
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SECTION  II 

AEROBALLISTIC  TESTING 


1.  Aeroballistic  Research  Facility 

The  Aeroballistic  Research  Facility'  (ARF)  is  an  enclosed  concrete  structure  used  to 
examine  the  exterior  ballistics  of  various  ffee-flight  configurations.  The  facility  contains  a  gun 
room,  control  room,  model  measurement  room,  blast  chamber,  and  an  instrumented  range. 
Figure  1  contains  an  illustration  of  the  ARF  and  an  interior  view.  The  range  atmosphere  is 
controlled  and  closely  monitored. 

The  200-meter  range  has  a  16-square  meter  cross-section  for  the  first  70-meters  and  a 
25-square  meter  cross-section  for  the  remaining  length.  The  range  has  131  locations  available 
as  instrumentation  sites  and  each  location  is  separated  by  1.5-meters.  Presently,  50  of  these 
sites  are  used  to  house  fully  instrumented  orthogonal  shadowgraph  stations.  At  each  of  these 
stations,  the  maximum  shadowgraph  window  (an  imaginary  circle  in  which  a  projectile  will 
cast  a  shadow  on  two  orthogonal  reflective  screens)  is  2-meters  in  diameter.  The  orthogonal 
photographs  of  the  model's  shadow  are  then  used  to  determine  the  spatial  position  and  angular 
orientation  of  the  test  model  at  each  of  the  50-instrumented  sites.  The  range  is  an  atmospheric 
test  facility  where  the  temperature  and  the  relative  humidity  are  controlled  to  22  ±1  °C  and  less 
than  55%  respectively. 

The  Comprehensive  Automated  Data  Reduction  and  Analysis  System  (CADRA)^  is 
used  to  read  the  film  and  calculate  the  trajectory.  The  film  is  digitized  using  a  high  resolution 
scanning  process.  Automated  image  processing  is  done  using  CADRA. 

A  chronograph  system  provides  the  flight  times  for  the  projectile  at  each  station.  These 
times  together  with  the  spatial  position  and  orientation  obtained  from  the  orthogonal 
photographs  processed  by  the  CADRA  system  provide  the  basic  trajectory  data  for  the 
subsequent  analysis.  These  discrete  times,  positions,  and  orientations  are  then  used  by  the 
Aeroballistic  Research  Facility  Data  Analysis  System  (ARFDAS)^  to  determine  the 
aerodynamic  coefficients  and  stability  derivatives  acting  on  the  model  during  the  observed 
flight. 
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a)  exterior  view  b)  interior  view 

Figure  1.  USAF  Aeroballistic  Research  Facility  (ARF),  Eglin  AFB,  FL 


2.  Models  and  Test  Conditions 

A  total  of  seven  configurations  were  designed  and  tested  for  this  effort.  All 
configurations  are  designed  to  have  an  equivalent  cross  sectional  area  (i.e.  reference  area).  In 
addition,  the  reference  length  used  in  the  data  analysis  was  common  for  all  configurations.  For 
missile  configurations,  the  reference  length  is  typically  the  diameter  of  the  projectile.  For  the 
non-circular  configurations,  an  “equivalent”  diameter  is  used  for  the  reference  length  such  that 
the  area  of  a  circle  based  upon  this  diameter  is  equivalent  to  the  non-circular  area.  However, 
since  the  cross-sectional  areas  were  equivalent  for  all  configurations,  the  reference  length  is  the 
same  for  all  configurations  as  well.  This  allows  for  one-to-one  comparison  of  the  data  of  each 
configuration  to  isolate  the  body  alone  influence. 

The  fin  planform  was  identical  for  all  configurations.  This  will  further  isolate  the 
effects  of  the  body  alone.  The  test  configurations  are  8.08  calibers  in  length  with  a  2.2  caliber 
ogive.  Figure  2  contains  a  schematic  illustrating  the  test  model  configurations  tested  during  this 
research  program.  The  model  dimensional  details  are  provided  in  the  following  sections  along 
with  the  aerodynamic  data.  A  fin  tab  was  attached  to  the  fm  tip  trailing  edge  of  one  fin  of  each 
model  in  order  to  determine  roll  orientation.  The  flight  data  covered  a  Mach  number  range  of 
0.75  to  3.5.  Table  1  lists  the  nominal  physical  properties  for  the  configurations  tested. 


Square  Trianoular 


Blended  Ellipse 


Figure  2.  Model  Cross  Section  ConHgurations 


Table  1:  Model  Physical  Properties 


Configuration 

Circular 

4  Fin 

Circular 

3  Fin 

0.8 

Elliptical 

4  Fin 

0.6 

Elliptical 

4  Fin 

Blended 

Elliptical 

4  Fin 

Square 

4  Fin 

Triangular 

3  Fin 

Diameter  {d),  mm 

17 

17 

17 

17 

17 

17 

Hfll 

Length  (L),  mm 

137.5 

137.5 

137.5 

137.5 

137.5 

Mass,  g 

844 

835 

916 

912 

853 

36.0 

34.6 

■6KW 

42.3 

39.8 

38.4 

1230 

■MiKM 

1290 

1410 

C.G.  location  {Xcg), 
mm  from  nose 

mil 

58.6 

64.8 

63.5 

jj^m 

H 

Number  of  ARF  Trials 

15 

15 

17 

15 

11 

13 

18 
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3.  Aerodynamic  Parameter  Identification 

Extraction  of  the  aerodynamic  coefficients  and  stability  derivatives  is  the  primary  goal 
in  analyzing  the  trajectories  measured  in  the  ARF.  The  process  is  summarized  in  Figure  3  and 
is  accomplished  by  using  ARFDAS^  ARFDAS  incorporates  a  standard  linear  theory  analysis'*’^ 
and  a  six  degree-of-freedom  (6DOF)  numerical  integration  technique*^.  The  6DOF  routine  in 
ARFDAS  incorporates  the  Maximum  Likelihood  Method  (MLM)’  to  match  the  theoretical 
trajectory  to  the  experimentally  measured  trajectory.  The  MLM  is  an  iterative  procedure  that 
adjusts  the  aerodynamic  coefficients  to  maximize  a  likelihood  function.  The  use  of  this 
likelihood  function  eliminates  the  inherent  assumption  in  least  squares  theory  that  the 
magnitude  of  the  measurement  noise  must  be  consistent  between  dynamic  parameters 
(irrespective  of  units).  In  general,  the  aerodynamics  can  be  nonlinear  functions  of  the  angle  of 
attack,  Mach  number,  and  aerodynamic  roll  angle.  ARFDAS  also  has  the  ability  to  analyze 
nonaxisymmetric  projectiles  using  a  “body  fixed”  aerodynamic  model.  The  “fixed  plane” 
method  for  the  analysis  of  axisymmetric  projectiles  is  given  in  Appendix  B  and  the  “body 
fixed”  method  is  described  in  Appendix  C. 

ARFDAS  represents  a  complete  ballistic  range  data  reduction  system  capable  of 
analyzing  both  symmetric  and  asymmetric  bodies.  The  essential  steps  of  the  data  reduction 
system  are  to;  (a)  assemble  the  dynamic  range  data  (time,  position,  attitude),  physical 
properties,  and  atmospheric  conditions,  (b)  perform  linear  theory  analysis,  and  (c)  perform 
6DOF  analysis  for  final  aerodynamics.  These  steps  have  been  integrated  into  ARFDAS  to 
provide  the  test  engineer  with  a  convenient  and  efficient  means  of  interaction.  At  each  step  in 
the  analysis,  permanent  records  for  each  flight  are  maintained  such  that  subsequent  analysis 
with  data  modifications  are  much  faster. 

Each  model  fired  in  the  ARF  was  initially  analyzed  separately,  then  combined  in 
appropriate  groups  for  simultaneous  analysis  using  the  multiple  fit  capability.  This  provides  a 
common  set  of  aerodynamics  that  match  each  of  the  separately  measured  position-attitude-time 
profiles.  The  multiple  fit  approach  provides  a  more  complete  spectrum  of  angular  and 
translational  motion  than  would  be  available  from  any  one  trajectory  considered  separately. 
This  increases  the  probability  that  the  determined  coefficients  define  the  model's  aerodynamics 
over  the  entire  range  of  test  conditions. 


ARFDAS  -  Aemballistic  Research  Facility 
Data  J^iysis  System 


Aerodynamic  Forces  &  Moments  vs. 
Mach  No.  &  Angle  of  Attack  &  Roll  Angle 

C  Xo  ^  XaA  ^Yp  a  ’  ^Yp  «3‘ 

a  <  a3  >  «5  >  ^mq  >^mq  2  '  ^npa>^«pa3>^npaS'  - 

y  ^nr  » ^nr  2  >  ^na  » ya  3*^nya3 . 

^tpy^l  pal » S^y  ya  3 . 


Figure  3.  ARFDAS  Aerodynamic  Parameter  IdentiHcation  Process 
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SECTION  III 

RESULTS  AND  DISCUSSION 


Aerodynamic  force  and  moment  coefficients  and  stability  derivatives  have  been 
extracted  from  each  set  of  free  flight  motion  data  measured  in  the  ARF.  The 
methodology  includes  both  linear  theory  and  six  degree-of- freedom  (6DOF)  matching  of 
the  observed  motions  to  determine  the  aerodynamics.  The  parameter  identification 
methodology  provides  a  best  match  to  the  experimentally  measured  motion  by 
determining  the  aerodynamic  forces  and  moments  acting  on  the  flight  vehicle  resulting 
in  the  measured  motion. 

The  following  sub-sections  present  the  aerodynamic  data  determined  from  the 
ballistic  range  trials  for  each  class  of  configuration:  circular,  elliptic,  square,  and 
triangular.  Appendix  D  contains  the  tabular  data  of  the  physical  properties,  range 
conditions,  and  6DOF  aerodynamic  results  for  each  configuration.  The  6DOF  results 
contain  both  the  single  and  multiple  fit  analysis  data.  It  is  believed  that  the  multiple  fit 
analysis  represent  the  best  estimates  of  the  configurations  aerodynamics  since  the 
combined  trajectories  contain  more  data  over  a  wider  range  of  angles  of  attack.  The 
plots  of  aerodynamic  data  presented  in  the  following  sub-sections  contain  hollow  and 
solid  data  symbols.  A  hollow  symbol  represents  the  result  of  a  single  fit  analysis  and  a 
solid  symbol  represents  the  result  of  matching  multiple  flight  trajectories  to  a  common 
set  of  aerodynamics.  It  should  be  noted  that  the  angular  motion  amplitude  for  most 
flights  was  very  small.  Such  small  amplitude  angular  and  swerve  motion  diminish  the 
accuracy  of  the  resultant  aerodynamic  coefficients  since  its  effect  on  matching  the 
measured  motion  becomes  less  significant. 

1.  Circular  Cross  Section  Configurations 

Two  circular  cross-section  configurations,  a  three-fin  and  a  four-fin,  were  tested 
to  provide  a  baseline  for  the  other  non-circular  configurations.  The  four-fin 
configuration  provides  a  baseline  for  the  elliptic  and  square  cross-section  configurations 
and  the  three-fin  provide  a  baseline  for  the  triangular  configuration.  The  two  circular 
cross-section  configurations  had  identical  bodies  and  fin  planform  dimensions.  Figure  4 
contains  a  schematic  of  these  two  configurations. 

Figure  5  presents  the  zero  yaw  axial  force  coefficient  (Cxo )  versus  Mach  number 
as  determined  from  the  flight  data  for  both  the  three-fin  and  four-fin  circular  cross 
section  configurations.  Since  the  only  difference  in  the  three  and  four  fin  models  is  the 
number  of  fins,  the  difference  in  drag  seen  in  this  figure  is  a  direct  result  of  the 
additional  fin.  The  data  indicate  there  is  about  a  5%  difference  in  drag  over  the  Mach 
number  range  from  0.6  through  3.2 


Q 


Figure  4.  Three-  Fin  and  Four-Fin  Circular  Configurations 

Figure  6  shows  the  normal  force  coefficient  derivative  (Cva)  versus  Mach 
number  for  the  three-fin  and  four-fin  circular  configurations.  Here,  the  additional 
normal  force  resulting  from  the  four  versus  three  fins  is  quantified.  There  is  about  an 
18%  difference  in  normal  force  coefficient  derivative  at  Mach  0.8.  The  difference 
increases  to  about  30%  at  Mach  1 .4  and  then  decreases  to  about  5%  at  Mach  3. 

Figure  7  contains  the  pitch  moment  coefficient  derivative  (C,„a)  as  a  function  of 
Mach  number.  The  difference  in  static  stability  for  three  versus  four  fins  is  clearly  seen 
and  is  accurately  determined.  The  difference  in  pitch  moment  coefficient  derivative  for 
the  three-fin  and  four-fin  configurations  is  about  20%  over  the  entire  Mach  number 
range. 


The  variation  of  pitch  damping  moment  coefficient  {Cmq)  versus  Mach  number  is 
shown  in  Figure  8.  Due  to  the  small  amplitude  motion,  it  was  difficult  to  accurately 
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determine  the  pitch  damping  moment.  However,  Figure  8  indicates  that  there  was  not  a 
significant  difference  in  the  pitch  damping  with  three  versus  four  fins. 


Figure  5.  Zero  Yaw  Axial  Force  Coefficient  versus  Mach  Number  -  Circular  Cross 

Section  Configurations 


Figure  6.  Normal  Force  Coefficient  Derivative  versus  Mach  Number  -  Circular 

Cross  Section  Configuration 
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Damping  Moment  Coefficient  ^  Pitch  Moment  Coefficient  Derivative 


2.  Elliptic  Cross  Section  Configurations 

Three  elliptic  cross-section  configurations  were  tested  in  this  effort.  The  first 
configuration  had  an  elliptic  cross  section  with  a  0.8  eccentricity  (ratio  of  minor 
diameter  to  major  diameter)  and  the  second  had  an  elliptic  cross  section  with  a  0.6 
eccentricity.  The  third  configuration  was  a  “blended”  cross  section  that  had  a  0.6 
eccentricity  on  the  top  half  and  a  0.8  eccentricity  on  the  lower  half  All  three 
configurations  contained  four  fins  and  are  shown  in  Figure  9.  As  noted  in  Table  1,  the 
center  of  gravity  for  each  configuration  was  different. 

Analysis  of  the  flight  data  was  done  using  the  6DOF  “body  fixed”  equations  of 
motion  within  ARFDAS^  (see  Appendix  C).  This  allowed  solving  for  unique  force  and 
moment  coefficients  in  the  pitch  (alpha)  and  yaw  (beta)  planes  (e.g.  Cma  ).  The 
four-fin  circular  cross  section  configuration  serves  as  a  baseline  for  the  elliptic 
configurations.  Since  all  configurations  tested  in  this  effort  have  a  common  cross 
sectional  area,  the  reference  length  and  area  are  common  to  all  configurations  and  this 
allows  for  direct  comparison  of  the  resultant  data.  However,  there  was  variation  in  the 
CG  locations  of  the  various  configurations  and  this  will  be  important  for  moment 
coefficient  comparisons. 

The  zero  yaw  axial  force  coefficient  {Cxo  )  versus  Mach  number  as  determined 
from  the  three  elliptic  cross  section  configuration  flight  data  is  contained  in  Figure  10. 
Also  included  for  comparison  are  the  4-fin  circular  cross  section  results.  Here  again,  the 
solid  symbols  represent  “multiple  fit”  data  versus  single  shot  analysis.  Although  there  is 
some  scatter  in  the  data  in  the  Mach  0.8  regime,  overall  the  data  appear  consistent  for  all 
configurations. 

Figure  11  shows  the  pitch-plane  (or  alpha-plane)  normal  force  coefficient 
derivative  (Cza )  versus  Mach  number.  The  data  of  Figure  1 1  show  the  increase  in  body 
lift  generated  with  the  elliptic  body  in  the  plane  of  the  major  diameter. 

Figure  12  shows  the  yaw-plane  (or  beta-plane)  yaw  force  coefficient  derivative 
(Cjp  )  versus  Mach  number.  The  yaw-plane  force  coefficient  derivative  results  have 
more  scatter  than  the  pitch-plane  force  coefficient  derivative  as  seen  in  Figure  1 1  due  to 
low  amplitude  motion.  However,  the  trend  shows  decreased  lift  in  the  yaw-plane 
(minor  diameter)  versus  the  pitch-plane. 
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Figure  9.  Elliptic  Cross  Section  Conflgurations 


Mach  Number 

Figure  10.  Zero  Yaw  Axial  Force  Coefficient  versus  Mach  Number  -  Elliptic  Cross 

Section  Configurations 


Figure  13  contains  the  pitch  (alpha-plane)  moment  coefficient  derivative  (C^a) 
as  a  function  of  Mach  number  as  extracted  from  the  ARF  flight  data  measuremOTts.  As 
indicated,  the  moment  references  were  the  individual  CG  locations  for  each 
configuration.  The  Mach  number  trends  are  similar  as  would  be  expected.  However,  to 
make  direct  comparisons  between  the  configurations,  the  moment  reference  should  be 
common  for  each  configuration.  Figure  14  contains  the  same  results  with  the  moments 
adjusted  to  a  common  body  reference  location  coincident  with  the  CG  location  for  the 
four-fin  circular  cross  section  model.  The  moment  adjustments  were  done  using  the 
following  equation^: 

C«a(58.6ww  Ref)  =  C„„(CG  Ref) +  ^ . ) 


Comparing  the  pitching  moment  results  for  all  three  elliptic  configurations  and 
the  circular  cross  section  (Figure  14),  the  differences  are  small  with  the  possible 
exception  of  the  transonic  Mach  regime.  This  would  imply  that  the  increase  in  Normal 
force  for  the  elliptic  shapes  was  coupled  with  a  forward  movement  of  the  center  of 
pressure  location  resulting  in  an  equivalent  pitching  moment. 


15 


Pitch-Plane  Force  Coefficient  Derivative  (C 


Mach  Number 


Figure  11.  Pitch-Plane  (Alpha-Plane)  Force  Coefficient  Derivative  versus  Mach 
Number  -  Elliptic  Cross  Section  Configurations 


Mach  Number 

Figure  12.  Yaw-Plane  (Beta-Plane)  Normal  Force  Coefficient  Derivative  versus 
Mach  Number  -  Elliptic  Cross  Section  Configurations 
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Figure  13.  Pitch-Plane  (Alpha-Plane)  Moment  Coefficient  Derivative  versus  Mach 
Number  -  Elliptic  Cross  Section  Configurations 


Figure  14.  Pitch-Plane  (Alpha-Plane)  Moment  Coefficient  Derivative  versus  Mach 
Number  with  Common  Moment  Reference  -  Elliptic  Cross  Section  Conflgurations 
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Figure  15  contains  the  yaw  (beta-plane)  moment  coefficient  derivative  (C„p)  as  a 
function  of  Mach  number  as  extracted  from  the  free-flight  data  measurements.  Here 
again,  as  in  the  case  for  the  pitch  plane,  the  moment  references  were  the  individual  CG 
locations  for  each  configuration.  The  yaw  moment  coefficient  derivatives  were  adjusted 
to  a  common  body  axial  location  and  the  results  are  plotted  in  Figure  16. 

Comparing  the  yawing  moment  results  for  all  three  elliptic  configurations  and 
the  circular  cross  section  (Figure  16),  the  differences  are  small  with  the  possible 
exception  of  the  transonic  Mach  regime  as  is  the  case  for  the  moment  in  the  alpha-plane. 
However,  the  yaw  force  (beta  plane)  showed  a  slight  decrease  as  compared  to  the 
circular  cross  section.  This  would  imply  that  the  decrease  in  yaw  force  for  the  elliptic 
shapes  was  coupled  with  a  rearward  movement  of  the  center  of  pressure  location 
resulting  in  an  equivalent  yawing  moment. 

The  damping  moment  coefficients  in  the  alpha  plane  {C^q)  and  the  beta  plane 
(C„r)  extracted  from  the  flight  data  are  plotted  in  Figure  17  and  Figure  18.  The 
qualitative  trends  show  a  small  decrease  in  the  damping  moment  in  the  beta  plane 
(minor  diameter). 


Figure  15.  Yaw-Plane  (Beta-Plane)  Moment  Coefficient  Derivative  versus  Mach 
Number  -  Elliptic  Cross  Section  Configurations 
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Pitch-Plane  Damping  Moment  Coefficient 


Mach  Number 


Figure  16.  Yaw-Plane  (Beta-Plane)  Moment  Coefficient  Derivative  versus  Mach 
Number  with  Common  Moment  Reference  -  Elliptic  Cross  Section  Configurations 


0.4  0.8  1.2  1.6  2.0  2.4  2.8  3.2  3.6 


Mach  Number 

Figure  17.  Pitch-Plane  (Alpha-Plane)  Damping  Moment  Coefficient  versus  Mach 
Number  -  Elliptic  Cross  Section  Configurations 


19 


Figure  18.  Yaw-Plane  (Beta-Plane)  Damping  Moment  Coefficient  versus  Mach 
Number  -  Elliptic  Cross  Section  Configurations 


The  center  of  pressure  can  be  found  via  the  following  equation^: 

^c.p.  _  ^CG  ^ma  ^  (2) 

L  L  C,,L 


The  center  of  pressure  for  the  pitch-plane  (alpha-plane)  was  computed  based  on 
the  resulting  normal  force  and  pitch  moment  coefficients  in  the  pitch-plane.  These 
results  are  presented  in  Figure  19  on  a  scale  consistent  with  the  total  body  length.  The 
center  of  pressure  for  the  yaw-plane  (beta-plane)  was  computed  in  a  similar  fashion 
based  on  the  resulting  yaw  force  and  yaw  moment  coefficients  in  the  beta  plane.  These 
results  are  presented  in  Figure  20  on  a  scale  consistent  with  the  total  body  length. 

In  Figure  21,  these  center  of  pressure  results  in  the  pitch-plane,  for  the  multiple 
fit  results  only,  are  plotted  on  an  expanded  scale.  This  provides  a  detailed  comparison 
of  the  most  accurate  results  showing  a  slight  forward  movement  of  the  normal  force 
center  of  pressure  (alpha  plane)  with  the  elliptic  configurations.  Figure  22  shows  the 
center  of  pressure  results  for  the  yaw-plane  on  the  same  expanded  scale  for  that  of  the 
pitch-plane.  In  general,  the  center  of  pressure  in  the  yaw-plane  is  about  5%  of  the  body 
length  aft  of  the  center  of  pressure  in  the  pitch-plane.  However,  the  data  is  not  refined 
to  make  an  accurate  quantitative  measurement. 
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Yaw-Plane  Center  of  Pressure  (%  body  length)  ^  Pitch-Plane  Center  of  Pressure  (%  body  length) 
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^'igure  19.  Pitch-Plane  (Alpha-Plane)  Center  of  Pressure  Location  versus  Mach 
Number  -  Elliptic  Cross  Section  Configurations 


A 

A 


CG  location  (4  fin  Circular  models}  / 

CG  location  (0.8  Elliptical  models)  ' 

CG  location  (Blended  Elliptical  models) 


CG  location  (O' 6  Elliptical  models) 


•  Circular  -  4  Fin 

A  0.8  Elliptical  •  4  Fin 

♦  0.6  Elliptical -4  Fin 

■  Blended  Elliptical  -  4  Fin 


Mach  Number 


Figure  20.  Yaw-Plane  (Beta-Plane)  Center  of  Pressure  Location  versus  Mach 
Number  -  Elliptic  Cross  Section  Configurations 
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Yaw-Plane  Center  of  Pressure  (%  body  length)  Pitch-Plane  Center  of  Pressure  (%  body  length) 


i’igure  21.  Pitch-Plane  (Alpha-Plane)  Center  of  Pressure  Location  versus  Mach 
Number  (enhanced  scale)  -  Elliptic  Cross  Section  ConOgurations 


Figure  22.  Yaw-Plane  (Beta-Plane)  Center  of  Pressure  Location  versus  Mach 
Number  (enhanced  scale)  -  Elliptic  Cross  Section  Configurations 
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3.  Square  Cross  Section  Configurations 

One  square  cross  section  configuration  was  tested  in  this  effort.  The 
aerodynamic  data  from  this  configuration  will  be  compared  with  the  four-fin  circular 
configuration  data.  Figure  23  depicts  the  four-fm  square  cross  section  configuration 
tested  in  this  effort.  It  should  be  noted  from  Figure  23  that  the  square  section  blends 
with  a  segment  of  the  ogive.  In  addition,  the  comers  of  the  square  are  rounded  as  shown 
in  Figure  23.  This  data  was  analyzed  using  the  fixed  plane  equations  of  motion  in 
ARFDAS. 


Figure  23.  Square  Body  Cross  Section  Configuration 


Figure  24  contains  the  zero  yaw  axial  force  coefficient  results  extracted  from  the 
flight  data.  The  supersonic  drag  for  the  square  cross  section  is  about  7  percent  higher 
than  the  circular  cross  section  but  at  subsonic  velocities,  there  is  very  little  difference  in 
the  data. 

Figure  25  contains  the  normal  force  coefficient  derivative  results  for  the  square 
cross  section  configuration.  Relative  to  the  four-fin  circular  cross  section,  the  normal 
force  is  larger  in  the  supersonic  regime,  only  slightly  larger  in  the  transonic  regime,  and 
about  the  same  subsonic. 

The  pitch  moment  coefficient  derivative  results  are  presented  in  Figure  26.  The 
square  cross  section  models  resulted  in  about  a  10%  increase  in  pitch  moment 
coefficient  derivative  versus  the  four-fin  circular  configuration  in  the  supersonic  regime. 
The  difference  is  smaller,  about  5%  in  the  subsonic  regime. 

The  pitch  damping  moment  coefficient  derivative  results  are  presented  in  Figure 

27.  Here  only  a  slight  change  from  the  four- fin  circular  configuration  is  noticed. 

The  normal  force  center  of  pressure  was  computed  based  on  the  extracted 
normal  force  and  pitching  moment  (see  equation  2).  The  results  are  plotted  in  Figure 

28.  Differences  in  center  of  pressure  location  between  the  two  configurations  are  small. 
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Coefficient  Derivative  (C/v„)  S  Zero  Yaw  Axial  Force  Coefficient  (C 


gure  24.  Zero  Yaw  Axial  Force  Coefflcient  versus  Mach  Number  -  Square  Cross 

Section  Configuration 


Figure  25.  Normal  Force  Coefficient  Derivative  versus  Mach  Number  -  Square 

Cross  Section  Configuration 


Figure  26.  Pitch  Moment  Coefficient  Derivative  versus  Mach  Number  -  Square 

Cross  Section  Configuration 


Figure  27.  Pitch  Moment  Damping  Coefficient  versus  Mach  Number  -  Square 

Cross  Section  Configuration 


25 


80 


40  -1- 
0.2 


♦ 


•  4  Fin  -  Circular 

♦  4  Fin -Square 


0-6  1.0  1.4  1.8  2.2  2.6  3.0  3.4 

Mach  Number 


Figure  28,  Center  of  Pressure  Location  versus  Mach  Number  (enhanced  scale)  - 

Square  Cross  Section  Conflguration 


4.  Triangular  Cross  Section  Configurations 

Figure  29  depicts  the  three-fin  triangular  cross  section  configuration  tested  in 
this  effort.  The  results  presented  are  compared  against  the  three-fm  circular  cross 
section.  It  should  be  noted  from  Figure  29  that  the  triangular  section  hlends  with  a 
segment  of  the  ogive.  In  addition,  the  comers  of  the  triangular  cross  section  are 
rounded. 


Triangular  Configuration 


Figure  29.  Triangular  Body  Cross  Section  Configuration 
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Initially,  the  extraction  of  the  aerodynamics  was  performed  with  the  body  fixed 
equations  of  motion  and  asymmetric  aerodynamic  model.  Analysis  of  the  data  was  also 
done  using  the  fixed  plane  model.  The  quality  of  the  match  to  the  measured  motion  was 
equivalent  for  both.  The  conclusion  is  that  the  differences  (e.g.  Cma  versus  C„p)  do  not 
have  a  measurable  effect  on  the  low  amplitude  motion  measurements  from  these  flights 
in  the  ARF.  Therefore,  additional  flights  with  induced  motion  for  larger  amplitudes  are 
required  to  adequately  define  the  aerodynamic  asymmetries. 

Figure  30  contains  the  axial  force  coefficient  results  extracted  from  the  flight 
rlata  The  transonic  and  low  supersonic  drag  for  the  triangular  cross  section  is  higher 
than  the  circular  cross  section. 

Figure  31  contains  the  normal  force  coefficient  derivative  results.  The  normal 
force  experimental  results  for  the  triangular  configuration  show  more  scatter  than  the 
circular  configuration  as  a  result  of  the  low  motion  amplitude.  However,  the  increase  in 
normal  force  resulting  from  the  triangular  cross  section  is  quantified. 

The  pitching  moment  coefficient  derivative  results  are  presented  in  Figure  32. 

The  pitching  moment  for  the  triangular  cross  section  is  larger  throughout  the 
Mach  range,  but  particularly  below  Mach  2.5. 

The  pitch  damping  results  are  presented  in  Figure  33.  Differences  between  the 
pitch  damping  coefficients  are  considered  within  the  measurement  noise  relative  to  the 
effect  of  matching  the  observed  motion. 

The  normal  force  center  of  pressure  was  computed  based  on  the  extracted 
normal  force  and  pitching  moment  via  equation  2.  The  results  are  plotted  in  Figure  34. 
The  center  of  pressure  for  the  triangular  cross  section  is  about  0.5  calibers  forward  as 
compared  to  the  circular  cross  section. 
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Normal  Force  Coefficient  Derivative  {Cfia)  Force  Coefficient  (C^ 
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igure  30.  Zero  Yaw  Axial  Force  Coefficient  versus  Mach  Number  -  Triangular 

Cross  Section  Configuration 
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Figure  31.  Normal  Force  Coefficient  Derivative  versus  Mach  Number  -  Triangular 

Cross  Section  Configuration 


Pitch  Damping  Moment  Coefficient  (C^)  ofj’  pitch  Moment  Coefficient  Derivative  (Q 


Mach  Number 


32.  Pitch  Moment  Coefficient  Derivative  versus  Mach  Number  -  Triangular 
Cross  Section  Configuration 


Mach  Number 


Figure  33.  Pitch  Damping  Moment  Coefficient  versus  Mach  Number  -  Triangular 

Cross  Section  Configuration 


80 


£ 

g  70 

s 

•D 

O 

X2 


3 

m 

01 

01 

%m 

(L 

O 

5 

s 

o 


60 


50 


▲ 

•  • 


A 


A 


A 


•  3  Fin  <  Circular 
A3  Fin  -  Triangular 


40  4- - ^ ^ ^ ^ _ i 

0.4  0.8  1.2  1.6  2.0  2.4  2.8  3.2  3.6 

Mach  Number 


Figure  34.  Center  of  Pressure  Location  versus  Mach  Number  (enhanced  scale) 
Triangular  Cross  Section  Configuration 
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SECTION  IV 

CONCLUSIONS 


This  report  documents  experimental  aerodynamic  test  results  of  a  variety  of 
cross  sectional  shapes;  all  of  which  have  the  same  cross  sectional  area.  Aerodynamic 
force  and  moment  coefficients  have  been  extracted  from  free-flight  data  measured  in  the 
AFRL  Aeroballistic  Research  Facility  (ARF).  The  range  of  Mach  numbers  covered 
during  these  trials  ranged  from  Mach  0.75  to  3.5  and  comprised  a  total  of  104  flights.  A 
total  of  seven  configurations  were  tested  and  evaluated:  four-fin  circular,  three-fin 
circular,  0.8  eccentric  elliptical  (four-fin),  0.6  eccentric  elliptical  (four-fin),  0.8/0.6 
eccentric  blended  elliptical  (four-fin),  four-fin  square,  and  three-fin  triangular. 

These  experimental  test  results  provide  a  good  comparison  of  the  shape  effects 
on  the  basic  aerodynamics  with  respect  to  axial  force,  normal  force,  and  pitching 
moment.  There  is  a  high  degree  of  confidence  in  the  aerodynamic  data  derived  from 
these  ballistics  range  tests.  There  was  repeatability  and  a  common  set  of  aerodynamics 
was  derived  from  multiple  flights. 

The  circular  cross  section  configurations,  a  4-fin  and  a  3-fin  configuration, 
established  a  baseline  reference  of  experimental  aerodynamics  for  comparison  to  the 
non-axisymmetric  configurations. 

With  respect  to  axial  force,  the  elliptical  cross  sections  did  not  exhibit  any 
increase  in  drag,  but  did  result  in  increased  normal  force.  However,  both  the  square  and 
triangular  cross  section  configurations  incurred  a  drag  penalty  to  achieve  an  increase  in 
normal  force. 

This  aeroballistic  test  program  has  provided  an  experimental  aerodynamic 
database  that  can  be  used  to  improve  aeroprediction  methodologies. 

Additional  testing  at  higher  angles  of  attack,  especially  for  the  elliptical 
configurations,  is  recommended.  The  test  flight  data  reported  here  do  not  contain  large 
enough  angular  motion  amplitude  to  accurately  quantify  aerodynamic  asymmetries. 
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APPENDIX  A  -  NOMENCLATURE 


6DOF 

Six  degree  of  freedom 

AFB 

Air  Force  Base 

AFRL 

Air  Force  Research  Laboratory 

AFRS 

Air  Force  Research  Shape 

ARF 

Aeroballistic  Research  Facility 

ARFDAS 

Aerohallistic  Research  Facility  Data  Analysis 
System 

C 

Celsius 

CADRA 

Comprehensive  Aerodynamic  Data  Reduction 
Analysis 

CG 

Center  of  gravity 

Ix 

Moment  of  inertia  about  x-axis 

Jxy 

Cross  product  moment  of  inertia 

fy 

Moment  of  inertia  about  y-axis 

Iz 

Moment  of  inertia  about  z-axis 

M 

Mach  number 

mbar 

milibar 

MLM 

Maximum  Likelihood  Method 

UF/GERC 

University  of  Florida  Graduate  Engineering  and 
Research  Center 

XcG 

Center  of  gravity  location  w.r.t  nose 

Xcp 

Center  of  pressure  location  w.r.t  nose 
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APPENDIX  B  -  FIXED  PLANE  AERODYNAMIC  MODEL 


1.  6DOF  -  Methodology 

The  aerodynamic  data  presented  in  this  report  obtained  using  the  “fixed  plane” 
6DOF  analysis  is  detailed  in  this  appendix.  Here,  the  equations  of  motion  are  derived 
with  respect  to  a  fixed  plane  coordinate  system.  The  x-axis  points  downrange,  the  y-axis 
points  to  the  left  looking  downrange,  and  the  z-axis  points  up. 

FIXED  PLANE  EQUATIONS  OF  MOTION 

u  =  gsinO-qw+rv-a  +F  Im 

CVi  X 

V  = -ru  -  rwianB  -  a  +F  Im 
cv  y 

iv  =  -?cos^  +  rvtan^  +  oM-a  +F  Im 

®  cw  z 

p^lU 

q  =  -r'^  Xaxi6 -{I ^1 1 y)rp  +  ml I y 

r  =  ^rtane  +  (/^ !  I  ^qp^n!  I  ^ 

Where  Ocm.  «cv»  ^cw  coriolis  accelerations  dependent  on  the  latitude  Xr 
and  azimuth  dg  of  the  range  and  rotational  rate  of  the  earth  (o^. 

BODY  FIXED  CORIOLIS  ACCELERATIONS 

a„  =  -2o)g  (y  sin  +z  cos  A,g  sin  ) 
ctg,  =  +2co^ (i sin - z cos Hg  cosSg) 

+2tUg  (jc  cos  Xg  sin  Sg+y  cos  Xg  cos  5g) 

Ucu  =  +c!cx  0  cos  cos  Osinif/  -  sin  0 

a^^-a^^sixi\i/-¥a^cosij/ 

=  +<^cx  ^ cos ^  sin esmy/  +  cos0 

Once  the  aerodynamic  forces  and  moments  (i.e.,  F^,  Fy,  F^,  I,  m,  n)  are  determined, 
the  solution  of  the  Equations  of  Motion  will  define  the  6DOF  flight  motion  with  respect 
to  the  fixed  plane  coordinate  system.  Since  the  position-attitude  measurements,  as 
acquired  from  the  ballistic  spark  range,  are  relative  to  the  Earth-fixed  coordinate  system, 
additional  transformation  equations  are  required.  These  transformation  equations  are 
shown  below  in  terms  of  the  fixed  plane  Euler  angles  (6,  v(;)  and  the  angle  of  rotation 
about  the  missile  axis  ((|»). 
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EARTH  FIXED  TRANSFORMATION  EQUATIONS  (FIXED  PLANE) 


x  =  u cos 0 cos y/  -vsintf/  +  ws\x\.9 cosy/ 
y  =  MCOS0sin;;/'  +  vcosy/  +  wsin^siin«/' 
i  =  -M  sin  ^  +  wcos  6 
e=q 

y/  =  r  I  cos  0 
^  =  p  +  r  tan0 


The  Equations  of  Motion  and  the  Earth  Fixed  Transformation  Equations 
numerically  integrated  using  a  fourth-order  Runge-Kutta  scheme. 

2.  Aerodynamic  Forces  and  Moments. 

The  fixed  plane  aerodynamic  forces  and  moments  are  defined  as  follows: 


where: 


F,  =-qAC^ 

_  -  .r  77  V  pd  -p:;  W 

y  -  q  [  ^  -f  ^ 

z  *2  L  Na  y  2V  V  V 

-  (Ca,<5  )  sin  -  (C^^0g )  sin  <j>] 

'=?4^c;,+c,] 

m  =  qAd[C„  -  +  ^C,„+£^  C.„  - 

n  L  ma  y  2y  J/' 

^aya  ^tia  ^mS^A  ^  “  ^mS^B 

-  jr  77  V  rd  p;  pd  p;  W 

H  —  Q^d\-~’C - h - C  “I - ^nnn  - 

'I  L  ma  y  2y  '”‘7  y 

—  W;  —  ^  — 

■*■  ^nya  ^na  ■*"  ^mS^A  ^  "*■  ^mS^B  COS(j) ] 


A  =  reference  area 
d  =  reference  length 
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q  =  dynamic  pressure  =  -^pV^ 


The  aerodynamic  coefficients  and  derivatives  are  assumed  to  be  nonlinear 
functions  of  Mach  number,  sine  of  the  total  angle  of  attack,  and  the  aerodynamic  roll 
angle.  This  assumption  is  made  in  a  general  sense  in  defining  a  generalized  aerodynamic 
math  model.  These  expansion  are  shown  as  follows: 

AERODYNAMIC  COEFFICIENT  EXPANSIONS  (FIXED  PLANE) 

Axial  Force  Coefficient 

C,  -  +  C,,/  +  C,,/  +  C^(M,  -M,)  +  oosN 

Normal  Force  Coefficient  Derivative 

C,a  =  f  ^  (M,  -  M, )  +  cos 

Magnus  Force  Coefficient  Derivative 

^Ypa  ~  ^Ypa  ^Ypa,^ 

Induced  Side  Force  Coefficient 

Sin  Nj 

Spin  Decay  Roll  Moment  Coefficient 
C»=C„+C»,.e^+C_(W,-M,) 

Slatic/lnduced  Roll  Moment  Coefllcient 

q  =C,^(J  +  C,^„2sinNV 

Pitching  Moment  Coefficient  Derivative 

=c„  +C.,y+C.^eUC.^(M,  -M,)+C,^,(M,  -M,f 
+  C„.  (CG  -  CG, )  +  £=  cos  Nr  +  C.„ 

Pitch  Damping  Moment  Coefficient 
=C„„  + 

Magnus  Moment  Coefficient  Derivative 

^/ipa  ~  ^npa  ^  ^"pccs  ^  ^npcun  i  ~  ^ r  ) 
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Induced  Side  Moment  Coefficient  Derivative 
=  C„ra  sin  sin 

Trim  Force  Coefficients 

^NS^A^^NS^B 

Trim  Moment  Coefficients 

^mS^A » B 

Out  of  Plane  Side  Moment  Due  to  Pitch 


The  aerodynamic  roll  angle,  y,  is  computed  by  transforming  the  fixed  plane 
missile  velocities  into  the  rolling  body  coordinate  system. 

Vj  =  vcos^  +  wsin^ 

=  -V  sin  ^  +  wcos  ^ 
y  =tan-'(v4/w4) 

The  sine  of  the  total  angle  of  attack  is  calculated  as  follows: 

£  = 


Slight  variations  in  the  center  of  gravity  (CG)  between  test  projectiles  (models) 
are  accounted  for  by  assigning  a  reference  CG  location  (CGj-)  and  making  an  appropriate 

correction  to  the  pitching  moment  coefficient  derivative.  The  pitching  moment 
coefficient  is  the  only  coefficient  of  which  slight  changes  in  CG  have  a  first  order  effect 
on  the  observed  motion. 

The  full  6DOF  equations  of  motion  portion  of  the  analysis  eliminates  the 
assumptions  of  Linear  Theory  by  retaining  all  cross  coupling  terms  and  allowing 
nonlinearities  both  as  functions  of  Mach  number  and  angle  of  attack.  In  addition,  the 
procedure  within  ARFDAS  allows  analysis  of  up  to  five  test  flights  simultaneously.  This 
provides  improved  accuracy  of  the  extracted  aerodynamics  and  their  nonlinearities  with 
angle  of  attack,  roll  angle,  and  Mach  number. 
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APPENDIX  C  -  BODY  FIXED  AERODYNAMIC  MODEL 


1. 6DOF  -  Methodology 

The  aerodynamic  data  presented  in  this  report  that  were  obtained  using  the  “body 
fixed”  6DOF  analysis  is  detailed  in  this  appendix.  Here,  the  equations  of  motion  are 
derived  with  respect  to  a  rotating  body  fixed  coordinate  system.  The  x-axis  points  down 
the  axis  of  the  body,  the  y-axis  points  out  the  left  side  of  the  body  looking  downrange, 
and  the  z-axis  points  up  with  respect  to  the  body.  The  body  fixed  coordinate  system  is 
rigidly  affixed  to  the  projectile  and  rotates  with  the  body  about  the  x-axis.  The  inertial 
frame  of  reference  is  the  earth.  It  is  assumed  the  earth  is  fixed  in  space  and  flat.  The 
body  fixed  equations  of  motion  are  given  as  follows  where  the  subscript  “h”  refers  to  the 
body  fixed  coordinate  system. 

BODY  FIXED  EQUATIONS  OF  MOTION 

=  g  Sine  -  q,,w^  +  +  — 

n 

h  =  Pb'^b  -  W?sin^cos^-a^,,ft  +-^ 

n  =  Qb^b  -  Pb^b  -  gcos<l>cos0-  a^f,  +  “ 

IJb  +  I^c^b  -  (4  +  h  -  h)hyPb^b  +  jJly  + 

Pb  ~  (II-  ) 

.  t  Vz,  +  (4  +  4  -  +  (h  (4  -  -  OPb^b 

^b  -  (/  7  -  7^  ) 

\^x^y  ^xy/ 

^xy^Pb  “  ”  ^y^Pb^b 

4=  ^ 

Where  Ucub,  dcvb.  and  OcWb  are  coriolis  accelerations  dependent  on  the  latitude  X/j 
and  azimuth  5g  of  the  range  and  rotational  rate  of  the  earth  ©g. 
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BODY  FIXED  CORIOLIS  ACCELERATIONS 


=  -2c0g(^siny?,yj  +zcos/l^  sin^^j) 
a^y  =  +2ft)g  (i  sin  Xr-z  cos  cosSr  ) 

=  +2cyg(;i:cosA^  sin^^j  +  ycos>l^  cos^^^j) 

^cub  -  ^cx  cos8cos\j/  +  a^y  cos^sin^ - sin0 
•  (sin  0  sin  ^  cos  -  cos  ^  sin  i//)  +  a ^y  (sin  0  sin  ^  sin  +  cos  ^  cos  ^j/') 

+  a^(sin^cos^) 

Ocwb  ~  ^cx  (sin  0  cos  ^  cos  ^  +  sin  ^  sinxj; )  +  (sin  0  cos  ^  sin  -  sin  ^  cos 

4-a^(cos^cos^) 

Once  the  aerodynamic  forces  and  moments  (i.e.,  Fx,  Fy,  F^,  I,  m,  n)  are  determined, 
the  solution  of  the  body  fixed  equations  of  motion  will  define  the  6DOF  flight  motion 
with  respect  to  the  body  fixed  coordinate  system.  Since  the  position-attitude 
measurements,  as  acquired  from  the  ballistic  spark  range,  are  relative  to  the  Earth-fixed 
coordinate  system,  additional  transformation  equations  are  required.  These 
transformation  equations  are  shown  below  in  terms  of  the  fixed  plane  Euler  angles  (9,  v|/) 
and  the  angle  of  rotation  about  the  missile  axis  (<})). 

EARTH  FIXED  TRANSFORMATION  EQUATIONS  (BODY  FIXED) 


The  Equations  of  Motion  and  the  Earth  Fixed  Transformation  Equations 
are  numerically  integrated  using  a  fourth-order  Runge-Kutta  scheme. 
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2.  Aerodynamic  Forces  and  Moments. 


where: 


F.,- 

Fyb- 


ttax  r  r  ^6  I  Pb^  c  -^1 

:  qAy  Cyo  V  ^  2V  V  ^  V  ^ 

Pb^  1 

2V  V  V^ 


V  --A\  r  r  Pb^  ( 

^zb~^^y~^Z<}~^Za  y  2y  ^^P°‘  V 

4-94^C,,  +  Q^+C,^] 

—  Ajy/-'  ,7=!  ^b  I  ^b^  7=!  .  Pb^  7^  Yb.i'^  ^b  i 

=  qAa[U^Q  +  y  ^  2V  2V  V 

-An  r  r  i^o-Z^r  4.Mr  3.  +  r  ^1 

=  qAai-Ly„Q  ”  '-np  ^  2F  F  F 


y4  =  reference  area 
d  =  reference  length 


nib- 
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=  dynamic  pressure  =  —  pF 


V  =V^b  +^6+^6 


The  aerodynamic  coefficients  and  derivatives  are  assumed  to  be  nonlinear 
functions  of  Mach  number,  sine  of  the  total  angle  of  attack,  and  the  aerodynamic  roll 
angle.  This  assumption  is  made  in  a  general  sense  in  defining  a  generalized  aerodynamic 
math  model.  These  expansion  for  the  body  fixed  equations  of  motion  are  shown  as 
follows: 


AERODYNAMIC  COEFFICIENT  EXPANSIONS  (BODY  FIXED) 
Axial  Force  Coefficient 

Cx  =  Cxo  +  Qa2  (yf  +  Cxfy  (yf  +  Cx.m  -  K)  +  Cxra2  ^OS  JV;. 

Normal  Force  Coefficient  Derivative 

Cza  =  Q«  +  Czas  (y)^  +  CNya^^ 
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Side  Force  Coefficient  Derivative 

Cyp  =  cosNy 

Magnus  Force  CoefHcient  Derivative 

^Ypa  ~  ^Ypa 

Induced  Side  Force  Coefficient 
^Yya  =  Cyya^s^  siniV;^ 

Spin  Decay  Roll  Moment  Coefficient 

Qp  =  Qp  +  Q/xx2^^  +  Qpm(Mi  -  M^) 

Static/Induced  Roll  Moment  Coefficient 

Pitching  Moment  Coefficient  Derivative 

?»a  =  C„ + C.^  +  Cz,(CG  -  CG,) + q„„s"  cos/V,- 


Yawing  Moment  Coefficient  Derivative 

C.H  =  C,f  +  C,^  +  C,^(CG  -  CG,)  +  c,„£"  cosAT^ 


Pitch  Damping  Moment  Coefficient 


c  -c  +  r 

^mq  ^mq  *  ^mqiX2  ^  y  ^ 


Yaw  Damping  Moment  Coefficient 


Qf  —  ^nr  +  Qm2  (  ^  ) 


Magnus  Moment  Coefficient  Derivative 


C  =C 

^npcL  ^npa 
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Induced  Side  Moment  Coefficient  Derivative 

Trim  Force  Coefficients 

^Z0>  ^YO 

Trim  Moment  Coefficients 

^mO’  Qo 

The  aerodynamic  roll  angle,  y,  is  computed  as  follows: 

y  =tan''(vjwj 

The  sine  of  the  total  angle  of  attack  is  calculated  as  follows: 


Iv^+w^ 

The  body  fixed  pitch  and  yaw  angles  are  defined  as  follows: 


Slight  variations  in  the  center  of  gravity  (CG)  between  test  projectiles  (models) 
are  accounted  for  by  assigning  a  reference  CG  location  (CGj-)  and  making  an  appropriate 

correction  to  the  pitching  moment  coefficient  derivative.  The  pitching  moment 
coefficient  is  the  only  coefficient  of  which  slight  changes  in  CG  have  a  first  order  effect 
on  the  observed  motion. 

The  full  6DOF  equations  of  motion  portion  of  the  analysis  eliminates  the 
assumptions  of  Linear  Theory  by  retaining  all  cross  coupling  terms  and  allowing 
nonlinearities  both  as  functions  of  Mach  number  and  angle  of  attack.  In  addition,  the 
procedure  within  ARFDAS  allows  analysis  of  up  to  five  test  flights  simultaneously.  This 
provides  improved  accuracy  of  the  extracted  aerodynamics  and  their  nonlinearities  with 
angle  of  attack,  roll  angle,  and  Mach  number. 
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Table  D-1:  4-Fiii  Circular  Model  Physical  Properties 
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Table  D-2:  4-Fm  Circular  Range  Conditions 
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Table  D-3:  4-Fin  Circular  6DOF  Aerodynamics  -  Single  Fits 
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Table  D-4:  4-Fin  Circular  6DOF  Aerodynamics  -  Multiple  Fits 
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2.  Three-Fin  Circular 
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Table  D-6:  3-Fm  Circular  Range  Conditions 
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Table  D-8:  3-Fin  Circular  6DOF  Aerodynamics  -  Multiple  Fits 
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3. 0.8  Elliptical 
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Table  D-10:  0.8  Elliptical  Range  Conditions 
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Table  D-11:  0.8  Elliptical  6DOF  Aerodynamics  -  Single  Fits 
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Table  D-12:  0.8  Elliptical  6DOF  Aerodynamics  -  Multiple  Fits 
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Table  D-13:  0.6  Elliptical  Model  Physical  Properties 
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Table  D-14:  0.6  Elliptical  Range  Conditions 
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Table  D-15:  0.6  Elliptical  6DOF  Aerodynamics  -  Single  Fits 
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Table  D-18:  Blended  Elliptical  Range  Conditions 
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Table  D-19:  Blended  Elliptical  6DOF  Aerodynamics  -  Single  Fits 
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Table  D-20:  Blended  Elliptical  6DOF  Aerodynamics  -  Multiple  Fits 
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Table  D-21:  Square  Model  Physical  Properties 
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Table  D-22:  Square  Range  Conditions 
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Table  D-23:  Square  6DOF  Aerodynamics  -  Single  Fits 
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Table  D>24:  Square  6DOF  Aerodynamics  -  Multiple  Fits 


RJ  0  ^ 
^  N 
CO  X  I 


s  s 

§IS 

u  u  u 


rvi  (N 
RJ  OJ 
(71  m  p4 

rH  ^  rH 

u  u  u 


to  ro  in 
(0  (d  n$ 
01  (T)  Di 


m  ro 
(d  Id  (d 
01  01  01 
N  0  0 


o  u  u 


g  §  § 

u  U  CJ 


(d  id  ro 
Oi  Id 

5h  d  tt 

u  u 

o 


VO  o  o 
H  ro  o 


OJ  H  O 
O  CM  O 


O  rH  O 
O  O  O 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

n 

o 

o 

o\ 

o 

o 

o 

o 

00 

o 

o 

VO 

o 

O 

CIV 

00 

d 

GO 

o 

o 

CO 

CM 

in 

w 

1 

ro 

• 

OJ 

CM 

1 

rH 

1 

r- 

o 

o 

o 

o 

o 

ro 

o 

o 

CM 

o 

o 

OV 

o 

O 

H 

1 

• 

ro 

« 

r- 

, 

* 

O 

, 

* 

VO 

00 

(N 

CN 

GO 

ov 

O 

o 

VO 

H 

1 

CM 

rH 

1 

rH 

1 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

O 

o 

O 

o 

o 

o 

• 

o 

o 

o 

o 

o 

o 

o 

n 

o 

o 

o 

o 

o 

O 

o 

• 

• 

(N 

• 

in 

* 

• 

o 

00 

d 

ov 

d 

d 

rH 

1^ 

o 

in 

OJ 

o 

ov 

in 

o 

VO 

ov 

o 

in 

o 

O 

in 

a\ 

o 

ov 

o 

ov 

CM 

VO 

ro 

» 

VO 

in 

ro 

• 

VO 

i/i 

rn 

n 

ov 

in 

in 

CN 

00 

ov 

VO 

r- 

oi 

ci 

ro 

rH 

ro 

VO 

n 

in 

ro 

rj 

00 

CM 

00 

ro 

o 

'-I* 

ro 

rH 

oj 

CM 

ro 

ov 

I> 

C'J 

o 

00 

VO 

r' 

r- 

CO 

'i* 

VO 

o 

o 

o 

rH 

in 

H 

H 

H 

rH 

o 

o 

ov 

o\ 

ov 

ov 

o 

o 

ov 

ov 

ov 

OV 

o 

o 

CO 

CO 

CO 

CO 

CO 

CO 

00 

in 

VO 

CM 

o 

in 

VO 

rH 

O' 

in 

in 

in 

o 

OJ 

o 

o 

vjf 

rH 

O 

rH 

o 

H 

o 

o 

OV 

o 

OV 

o 

ov 

o 

o 

cn 

o 

OV 

o 

ov 

o 

o 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

68 


Table  D-25:  Triangular  Model  Physical  Properties 
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Table  D-26:  Triangular  Range  Conditions 
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Table  D-27:  Triangular  6DOF  Aerodynamics  -  Single  Fits 
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Table  D-28:  Triangular  6DOF  Aerodynamics  -  Multiple  Fits 
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